Stability of a Fully Polarized Ultracold Fermi Gas near Zero-Crossing of a p-wave 

Feshbach Resonance 
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We consider a fully polarized ultracold Fermi gas interacting through a p- wave Feshbach resonance. 
Using a two-channel model, we find the effective potential at the point where the p-wave scattering 
length goes to zero. Here the effective interaction provides attraction and one can therefore ask 
about the stability of the system. We calculate the energy density of the system in the Thomas- 
Fermi approximation, determine the profile of the gas, and the critical number of particle in the 
system as function of the relevant interaction parameters. The instability can be inferred from a 
simple breathing mode argument which explains the scaling found numerically. The critical particle 
number turns out to be extremely large unless the external confinement is very tight. We therefore 
conclude that the effect is insignificant for standard trapping potentials and that the magnetic dipole 
interaction is the important term at zero scattering length. However, for tight confinement as in an 
optical lattice higher-order corrections can become important. 

PACS numbers: 03.75.Ss,67.85-d,37.10. Jk 



I. INTRODUCTION 

Feshbach resonances in ultracold atomic gases that al- 
low tuning of interactions are a continuing source of ex- 
citement. The use of s-wave resonances have led to great 
insight into f.x. the BCS-BEC crossover in the two- 
component Fermi gas Recently, also p-wave reso- 
nances have been explored in experiments with 40 K 041] 
and 6 Li and p-wave Feshbach molecules [|| have 

been created d, [1(3, ITTj |. This opens up the possibility 
of creating exotic p-wave superfluids for which several 
classical and quantum phase transitions have been pre- 
dicted theoretically |13 - [l9f . Simultaneously, there has 
been a lot of recent interest in using Feshbach resonances 
to tune the scattering length to zero [20l - |25j . These ex- 
periments work with Bose condensates and are able to 
probe residual magnetic dipole-dipole interactions at this 
zero-crossing (23|. However, as discussed in [26l - |28j there 
are also higher-order interaction terms present at zero- 
crossing that can influence the stability of the condensate 
for certain resonances. 

In this work we want to explore the effects of higher- 
order interaction terms on Fermi gases as the p-wave scat- 
tering length is tuned to zero. In particular, we study the 
fully polarized single-component system where s-wave in- 
teractions are Pauli suppressed and the p-wave term be- 
comes important. Such highly polarized Fer mi gase s are 
currently under intense study experimentally |29l - [34j . We 
will study some bulk properties of the trapped system 
through the local-density approximation. We note that 
the consequences of p-wave resonances have also recently 
been explored in traps for two- and three-particle systems 

As we will show below, the higher-order correction 
term due to the p-wave Feshbach resonance is attractive 
and can therefore serve to destabilize the system when 
the particle number becomes large. We consider the sta- 
bility toward collapse within mean-field theory and find 



that a necessary condition for the higher-order terms to 
be important is that the external confinement be very 
tight, i.e. the trapping length scale should be of order 
0.03 /jm or less. This is of course much smaller than typ- 
ical traps, however, it might be relevant for optical lattice 
setups. 

The paper is structured as follows. In Sec. [TT] we in- 
troduce and discuss a two-channel model for the p-wave 
resonance that is needed to describe the higher-order in- 
teraction term. We consider the full T-matrix at low 
momenta, find the effective potential, and then match 
the parameters to low-energy p-wave scattering. Then in 
Sec. IIIII we derive mean-field equations for the gas within 
the Thomas-Fermi or local-density approximation and 
discuss the density profiles. In Sec. lIII Al we present some 
analytical stability conditions and calculate the critical 
number of particles in the system as a function of the in- 
teraction strength. We compare and discuss experimen- 
tally relevant parameters in Sec. IIVI In the concluding 
Sec.|V]we summarize our findings and consider the rela- 
tive importance of higher-order terms from the resonance 
to the magnetic dipole-dipole interaction. Here we again 
find that very tight confinement must be used for the 
higher-order correction to be of comparable magnitude. 



II. p-WAVE EFFECTIVE POTENTIAL 



Since we are interested in higher-order interaction 
terms at zero-crossing we use a two-channel model to 
describe the p-wave Feshbach resonance. This model is 
obtained from quantum defect theory [38], [39| . If we as- 
sume that the energy of the bound state varies linearly 
with magnetic field, B, then the p-wave phase-shift sat- 
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where Bo is the resonance position, AB the width, Aji 
the difference in magnetic moment between the open and 
closed channel, and m the reduced mass of the scattered 
atoms. The function ftg(k) represents the background 
scattering from the long-range van der Waals (vdW) in- 
teraction. We will assume that E <C E v dw, where E 
is the scattering energy and E v dw is the energy scale of 
the vdW interaction. In this case we can ignore effec- 
tive range corrections from the background term and we 
have fbg(k) ~ ( a b g ) 3 i where a p g is the background p- wave 
scattering length. 

We can now write the on-shell open-open p-wave T- 
matrix for the two-channel model at low energy as 



T p (k) = 



4^ 2 « ) 3 
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At this point one could attempt an effective-range expan- 
sion and find the effective potential through the energy 
shift method [4(| • Unfortunately, at zero-crossing where 
AB = B — Bo the coefficients of such an expansion are all 
divergent as discussed for the s-wave case in [28]]. How- 
ever, since the full T-matrix is available we can follow 
the strategy of 28] . We go to zero-crossing and expand 
the T-matrix to lowest order in k to obtain 



Tp(fc) = rfK £ )Mo fc4 + o(fc6) 



(3) 



where we have introduced the background effective range 
r e0 = a~ar " ■ Notice that r P < in this zero-range 
model 44]. One can now solve the Lippmann-Schwinger 
equation to obtain the effective potential. Since we 
are interested in ultracold fermions we take the long- 
wavelength limit, i.e. we let the cut-off momentum in the 
solution go to zero 28] . In that case we find the immedi- 
ate answer V — T and we merely have to transform the 
momentum-space expression into coordinate space. In 
our case we need a second-order p-wave zero-range term 
as outlined in 40]. The properly symmetrized effective 
potential is 
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where g p 



(i) 



4 rf WJM 



and r is the relative distance 



of the two atoms. The superscript on the coupling in- 
dicates that we consider the first correction term to the 
scattering length approximation [401 ] . 



Since r p < we have 



< 0. As we show in the 
next section, the negative coupling leads to an attractive 
interaction term in the energy density. We can therefore 
ask about the stability of a single-component Fermi gas 
interacting through a p-wave Feshbach resonance around 
zero-crossing. This is determined by the interplay of ki- 
netic energy (or Fermi pressure) which tends to expand 
the gas and the attractive term above which causes con- 
traction. If the latter is dominant we expect the system 
to collapse. 

The collapse proceed through recombination events 
into deeply bound dimer states in the van der Waals po- 
tential between the atoms. In the case of a fully polar- 
ized Fermi gas, s-wave collisions are suppressed due to 
the Pauli principle. The leading order loss mechanism 
is therefore in the p-wave channel (3(| [U |42| . Close to 
a p-wave Feshbach resonance one can estimate the de- 
cay rate of the gas on both the BCS (a p < 0) and BEC 
(cip > 0) sides through the background parameters of 
the interaction as [42| T ~ ^mm ; where R v aw is the 
van der Waals length of the atoms [43| and I = ?i -1 / 3 
is the interatomic distance when n is the density of the 
gas. Based on this estimates one expects the lifetime of 
the polarized gases with a p-wave resonance to be much 
shorter than that of two-component gases with an s-wave 
resonance 42]. In the current work we do not consider 
the immediate vicinity of the resonance but rather zero- 
crossing where a p = 0. However, we expect the same 
estimates to hold in this case as well. As we will see be- 
low, the collapse near zero-crossing is associated with an 
instability of the breathing mode which entails a change 
in the radial extension of the cloud. We thus imagine 
the collapse taking place through a decrease in cloud size 
and thus in I with a corresponding increase in T. The 
instability should therefore manifest itself as a reduction 
of the lifetime of the system close to the critical number 
of particles calculated below. 



III. THOMAS-FERMI ENERGY DENSITY 

The energy density of a single-component homoge- 
neous Fermi gas with the effective interaction of Eq. ([4| 
can be written in terms of the Fermi momentum, kp. 
Here we use the Thomas-Fermi or local-density approx- 
imation to write the energy density in an external trap- 
ping potential in terms of the local Fermi momentum 
fcp(x) as 



h 2 k F (xf 
207T^m 



3^ 



— Ve X t(x)k F (xf + ^-—k F (x) W , 

Gn 2 ea; H 2800tt 4 to y ' ' 



(5) 



where the last term is obtained from the matrix ele- 
ments of Eq. in the Fermi gas (see and refer- 
ences therein for details). The local density is given by 
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FIG. 1: Spherically symmetric density profiles p(r) for 
p/pmax = 0.25, 0.5, 0.75, and 1.0 (from bottom to top) where 
fJ-max is the critical chemical potential calculated for the 40 K 
resonance discussed in the text and trap size b = 0.1pm. The 
solid (black) lines are obtained from Eq. ([8|, whereas the 
dashed (red) lines are the non-interacting densities for the 
same particle number. The particle numbers are N = 1.9- 10 3 , 
1.5 ■ 10 4 , 5.1 ■ 10 4 , and 1.2 • 10 5 (bottom to top). 



p(x) = kp(x) 3 /6ir 2 . The formalism we use is at zero tem- 
perature. This should be sufficient for the stability anal- 
ysis as experiments have temperatures of T/Tf < 0.05 
with Tp the Fermi temperature [HI, [34| . In the follow- 
ing we will assume that the external confinement has the 
form of an isotropic harmonic oscillator 
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with frequency u> and oscillator length 6, although we 
note that isotropy is not essential for the stability analysis 
presented below. 

We now determine the ground-state density of the 
single-component Fermi system by minimizing the energy 
for a given number of particles. To this end we introduce 
the chemical potential, p, and perform the variation 



d 3 x (£tf{x) — pp{x)) 
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at each point x. We then arrive at the equation 
=\{k F { X )bf 
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where explicit oscillator units have been introduced. For 
each x this represents a polynomial equation which can 
be solved for kp(x). The particle number, N, can be 
adjusted by varying p to ensure that N = J d 3 xp(x). 

In Fig. (p} we show some typical density profiles ob- 
tained by solving Eq. (|SJ) for p/p max — 0.25, 0.5, 0.75, 



and 1.0 (from bottom to top), where p m ax is the criti- 
cal chemical potential discussed in the next section. The 
full (black) lines are obtained from Eq. ([8]), whereas the 
dashed (red) lines are the corresponding non-interacting 
systems in the trap with the chemical potential adjusted 
to yield the same particle number. The interacting den- 
sities are in general higher in the center of the trap as 
expected from the attractive nature of the interaction 
at zero-crossing. In fact, close to the critical number of 
particles we see a clear difference in the density profiles, 
whereas for the lower particle numbers the interacting 
and non-interacting densities are very similar. The effects 
of higher-order interactions around zero-crossing can thus 
in principle be probed through the density profile of the 
gas. In the momentum distribution there should also be 
a signature of the squeezing of the cloud due to the at- 
traction. 



Critical Particle Number 



From the structure of Eq. {5} and the fact that Qp < 
we see that there will be a maximum Fermi momentum, 
kmax, and a corresponding p m ax beyond which no solu- 
tion can be found for a given x_(as discussed for the lowest 
order p-wave interaction in [40]). This means that there 
will be a critical number of particles, N max , so that for 
N > N max the system is unstable in mean-field theory. 
We now determine these quantities. 

In the center of the harmonic trap we have V ex t (0) = 0, 
and away from the center p/hui — (x/b) 2 /2 is monotoni- 
cally decreasing toward the edge at x/b = yj2p/hu>. The 
maximum Fermi momentum k max is therefore obtained 
by finding the turning point of the right-hand side of 
Eq. ©. This yields 
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where we have introduced a = (a% /b) \r^ /b\ to reduce 
the notation. In terms of density, the stability condition 
is a 3 ^ 5 p(x)b 3 < 0.07, which must hold everywhere in the 
trap for the system to remain stable. As discussed in 
40] , this condition is independent of the trap geometry 
and the use of an isotropic harmonic trap is therefore not 
essential for the stability analysis. 

From the values of k max above we can calculate the 
density profile for p = p max and determine N max by 
simple integration. In Fig. ^ we plot N max as a func- 
tion of the interaction strength parametrized by a on 
a log-log scale. We indeed see that the critical num- 
ber of particles grows rapidly with decreasing interaction 
strength. The (red) squares denote some estimates for a 
realistic p-wave resonance in 40 K (to be discussed below) 
with the upper one corresponding to b = lyum and the 
lower one having b = 0.1/im. In the former case we find 
N m ax ~ 10 11 whereas the tighter trapped system has 
Nmax ~ 10 5 which is in the regime of experiments with 
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FIG. 2: Critical particle number, N max , as function of the 
interaction strength parametrized by (a^ g /b) 4 \r^ /b\. The 
(black) circles indicate the decades. The (red) squares are 
for the 40 K resonance discussed in the text with oscillator 
length 6 = 1/im (left) and b = 0.1/im (right). 



-1.2 



polarized Fermi gases |29h34| . From the slope in Fig. ([2]) 
Decreasing the trap length by one 



oc a 



we find N„ 

order of magnitude is therefore seen to decrease N max by 
six orders of magnitude (in accord with the difference 
between the (red) squares in Fig. ([2])). 

We now proceed to show how the scaling of N max 
with a is related to the instability of the monopole or 
breathing mode of the system. To see this we employ the 
method of collective coordinates |43[ , although the same 
result can be found through the sum- rule approach [45| . 
We assume that the cloud maintains its shape but allow 
for variations of the spatial extent parametrized by the 
radius R. Considering small deviations from equilibrium 
we can derive a virial theorem of the form 



-2E K + 2E V - 7Ej = 0, 



(10) 



where the kinetic, Ek 1 harmonic trapping, Ey, and in- 
teraction, Ei, parts are the integrated equivalents of the 
terms in Eq. With this relation, we can relate the 
frequency of the breathing mode, ojb, to the energies 
through 



■ 35^ 
2 E v 



(11) 



The instability occurs since Ej oc —a < 0. For suffi- 
ciently large a the breathing mode frequency will there- 
fore become imaginary indicating an instability toward 
collapse when Ei/Ey < —8/35. 

The dependence of the instability condition on a and 
N can be found using the simple non-interacting density 
k F (x) = ^/2^/huj - x 2 /b 2 /b. We find E v = c v N 4 ^ 
and Ej = — ac/A^ 13 / 6 , where cy and cj are constants. 

5/6 

The instability condition is therefore aNmax — c, with 
c a constant of order one. We thus recover the scaling 



N m ax oc a as found numerically from the interacting 
densities above. 



IV. COMPARISON TO EXPERIMENTS 

We now estimate the parameters needed in our sce- 
nario for p-wave Feshbach resonances in 40 K 0]. These 
data have been fitted to the p-wave resonance model in- 
troduced above for all magnetic substates mi = 0, ±1 of 



|) hyperfine state in 37 1. The pa- 



the |/ 

rameters needed here are similar for all mi and for con- 
creteness we focus on the m; = state. The resonance is 
found at B = 198.85G with a width of AB = -20.342G, 
a background scattering length of a^ g — — 104.26ao (with 
ao the Bohr radius), and magnetic moment difference 
Afi = 0.067/is where /is is the Bohr magneton. For 
these values we find r^ = — 914.72ao From AB alone 
the resonance appears rather broad and a small value 
of the effective range would be expected. However, the 
small value of A/i compensates and increases r p eQ . In fact, 
p-wave resonances are in general very narrow as discussed 
in @. 

For the interaction parameter introduced above we 
have 



1.08 



io n (t: 



4.48 • 10" 



1 fim 



(12) 



and we thus immediately see that the tightness of the 
trapping is the decisive factor. Comparing to Fig. © we 
see that typical experimental trap lengths of b = 1/im will 
yield critical particle numbers that are several orders of 
magnitude above current experimental numbers (upper 
(red) square). However, as mentioned above, an order 
of magnitude decrease in the trap length will bring the 
number into the experimental range (lower (red) square). 
The density profiles for b — 0.1/im are shown in Fig. ([T]) 
together with the non-interacting densities and at the 
critical number we do see considerable deviations within 
the mean-field Thomas-Fermi picture. We conclude that 
tight trapping with b < 0.1/zm is necessary in order to 
experimentally observe any effect of higher-order inter- 
actions on the stability of the system around the zero- 
crossing of realistic p-wave Feshbach resonances. 

In this work we have exclusively employed the zero 
temperature formalism. One can of course extend the 
semi-classical approach to finite temperature to get a bet- 
ter description of experiments. Here we will be content 
with qualitative estimates of the effect of temperature 
which can be obtained by considering the non-interacting 
Fermi gas in a harmonic trap. As shown in the earlier 
work of [46j j , the intuitively clear effect of finite tempera- 
ture is to extend the density profile of the gas (see Fig. (3) 
of |46|). In turn, the local Fermi wavevector is reduced 
and likewise the effect of interaction in our Eq. We 
therefore expect the critical number to increase with tem- 
perature. However, as noted above experiments can cool 
to only 5-10% of the Fermi temperature and for this level 
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of degeneracy the deviation from the T = density pro- 
file should be very small ,46] . This justifies the use of the 
zero temperature formalism. 



V. CONCLUSIONS 

We have studied the effects of higher-order interactions 
in single-component Fermi gases interacting via a p-wave 
Feshbach resonance as the scattering length goes to zero. 
Using the low-energy limit of the full T-matrix we found a 
formula for the effective potential with strength propor- 
tional to the background parameters of the resonance. 
Implementing this effective interaction in a local-density 
Thomas- Fermi approximation for the Fermi gas in an iso- 
tropic harmonic oscillator trap, we found an equation for 
the density profile. Since the effective interaction around 
zero-crossing is attractive, the system has a critical par- 
ticle number beyond which it becomes unstable accord- 
ing to mean-field theory. We gave an analytic expres- 
sion for the stability condition and numerically calcu- 
lated the critical particle number for relevant experimen- 
tal parameters. The density profile was compared to the 
non-interacting profile and considerable deviations found 
near the critical particle number. In comparison to exper- 
imental parameters, we find that the external potential 
must provide a very tight confinement with trap length 
below 0.1pm for the critical number to be at the value of 
experiments. 

At zero-crossing there are other residual forces in the 
gas . In particular, as was recently found in Bose gases 
|22l |23| , the magnetic dipole-dipole interaction must be 
taken into account. We therefore have to consider the 
relative strength of this interaction to the higher-order 
p-wave term discussed in the current work. Comparing 
the couplings we find [43[ 



U (i) 

9 P 
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where the trap length scale b appears since the effective 
potential carries four derivatives. For the ratio to be 
of order one, a trapping length of b ~ 0.03pm is needed. 
With normal traps the magnetic dipole-dipole interaction 
is thus much stronger at zero-crossing. Curiously, we note 
that with an optical lattice of some 30 lattice sites in a 



beam waist of lpm the lattice site size is roughly 0.03pm. 
Exploring the effects of higher-order interactions in op- 
tical lattices is therefore a promising direction in which 
the current work can be extended. Of course, this also 
holds true for higher-order effects near zero-crossings for 
s-wave Feshbach resonances and work in this direction is 
also in progress. 

The attractive nature of the higher-order correction 
could perhaps also lead to p-wave superfluidity at zero- 
crossing. We can make an estimate of the critical tem- 
perature for this phase if we assume the same functional 
form as for s-wave superfluidity in two-component gases. 
As we expect it to increase with interaction strength and 
with Fermi wavevector, T c ~ Tp exp y— a ( fc ^ b )t, ^ ■ In the 

center of the trap, kp = k max , and from Eq. Q we 
find T c ~ T F exp(-3/107r) » 0.9T F which is a large 
value compared to experimentally realizable tempera- 
tures. However, factors of order 1-10 both inside and 
outside the exponential can severely reduce T c as seen in 
the s-wave case [43J, and we also have to keep in mind 
that Uf decreases toward the edge of the trap. A po- 
tentially much more severe hindrance is the fact that the 
estimates of the decay rate discussed above also suggest 
that the p-wave superfluid does not have time to form 
before it decays [421 ]. 

The above estimate demonstrates that the dipole- 
dipole interaction is overwhelmingly dominant in large 
traps, and an optical lattices is needed for the correc- 
tion to matter. Transition temperatures could be raised 
by manipulating the lattice symmetries as found in [47j 
and a stabilizing mechanism in the lattice was also re- 
cently proposed (48|. We aim to explore the influence 
of the higher-order corrections on such scenarios as well. 
We caution, however, that despite detailed experimental 
studies of p-wave Feshbach resonances [1, H| , an impor- 
tant parameter related to the Feshbach channel coupling 
cannot yet be deduced from the data and the possibility 
of producing p-wave superfluids is therefore still hard to 
estimate [1J]. In contrast the p-wave resonance model 
used above only requires background parameters and the 
magnetic moment difference. The evaluation of instabil- 
ity in mean-field theory is therefore a simpler problem. 
We still have to remember that our parameters where 
obtained from a fitting procedure of the data in Q to 
a two-channel model which is accurate to within a few 
percent [37| . 
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